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I. INTRODUCTION 



The languages of Hamiltonian and Lagrangian mechanics permeate much of fluid and 
plasma physics. On the one hand, Hamiltonian methods provide powerful analytical and 
numerical tools for the investigation of the motion of charged particles in complex electro- 
magnetic fields and the linear and nonlinear stability of plasma and fluid equilibria. On the 
other hand, Lagrangian (or variational) formulations of fluid and plasma physics provide 
clear guiding principles for the systematic dynamical reduction of fluid and plasma dynam- 
ics through the asymptotic elimination of fast degrees of freedom and the construction of 
adiabatic invariants either by Lagrangian averaging or by Lie-transform methods. More- 
over, variational formulations facilitate the derivation of important conservation laws (e.g., 
energy-momentum, angular momentum, and wave action) through the Noether method. 

The present paper summarizes the oral talks presented at the Mini- Conference on Hamil- 
tonian and Lagrangian methods in fluid and plasma physics held at the 2002 annual meeting 
of the Division of Plasma Physics of the American Physical Society. The paper is organized 
by topic and not the chronological order in which the talks were given. The organization 
is, therefore, as follows: Hamiltonian methods (Sec. II), Lie-transform methods (Sec. Ill), 
Lagrangian methods (Sec. IV), Lagrangian chaos (sec. V), and a Summary of the Mini- 
Conference (Sec. VI). Because this is a summary of a topical mini-conference, and not a 
review paper, no attempt has been made to provide a comprehensive list of references. In- 
stead, references to the relevant review papers have been provided, and citations to original 
sources can be found there. 

II. HAMILTONIAN METHODS 

Applications of Hamiltonian methods relie on the existence of a Hamiltonian function 
(or functional) if(£,t), which depends on a set of coordinates (or fields) ( k {k = 1, ...,N), 
and a Poisson bracket structure { , } such that Hamilton's equations are expressed as 
d( k /dt = {( k , H(£,t)}. Here, the Poisson bracket must satisfy (i) the antisymmetry prop- 
erty {/, g} = - {g, /}, (ii) the Leibnitz property {(f g), h} = f {g, h} + {/, h} g, and 
(iii) the Jacobi condition {/, {g, h}} = -{g, {h, /}} - {h, {/, g}}, where (f,g,h) 
are arbitrary functions (or functionals) . Whether the coordinates £ are canonical or not 
determines additional features of the Hamiltonian dynamics. Although the Poisson bracket 
structure of a canonical Hamiltonian formulation has a simple form, canonical coordinates 
suffer from the problem of being unphysical in most cases (e.g., they are gauge dependent). 
The exploitation of noncanonical Hamiltonian formulations was pioneered by the work of 
Littlejohn on the Hamiltonian theory of guiding-center motion [TJ. The physical trans- 
parency of noncanonical coordinates comes at the cost of the complexity of the Poisson 
bracket structure, although its derivation is solidly grounded mathematically. An impor- 
tant consequence of the noncanonical Poisson bracket structure is the existence of Casimir 
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invariants C(C), for which {C(Q, -F(C)} — for any function F(Q. For example, the non- 
canonical Lie-Poisson bracket structure { , }/ for the Vlasov equation is an antisymmetric 
bilinear operator on the space of Vlasov distributions f(z) in six-dimensional phase space 
defined as [2] 



where F[f] and G[f] are two arbitrary functionals and { , } is the Poisson bracket structure 
on six-dimensional phase space. Casimir invariants for the Vlasov equation are of the form 
C[f] = I x{f{ z )) d 6 z, where x(f) * s an arbitrary function of / and {C[/], F[f]}f = for all 
functionals F. Casimir invariants [3] play an important role in defining relative Hamiltonian 
equilibria and in investigating linear and nonlinear stability. The present Section focuses 
on some applications of Hamiltonian methods in fluid and plasma physics. 

Phil Morrison, the Mini- Conference keynote speaker, addressed the important question 
associated with the significance of the existence of a Hamiltonian structure or a variational 
(Lagrangian) structure for any dissipationless dynamical system ^j. First, the existence of 
a variational structure serves as a guiding principle for deriving exact and reduced dynami- 
cal equations j3] and their associated conservation laws. Several variational formulations for 
fluids and plasmas exist; the variational formulations for the Vlasov equation are reviewed 
by Ye and Morrison [Sj while the variational formulations for ideal fluids are reviewed by 
Holm, Marsden, and Ratiu [DJ. Next, from a practical point of view, the existence of a non- 
canonical Hamiltonian structure reveals the existence of Casimir invariants for a continuous 
Hamiltonian system and, thus, allow the investigation of the stability of relative Hamilto- 
nian equilibria [3]. Through the use of Hamiltonian methods, the stability properties of 
numerous dissipationless kinetic, hydrodynamic, and magnet ohydro dynamic systems have 
been investigated through the Lyapunov [7j, Energy-Casimir [S], and Dynamical Accessi- 
bility methods. Another significant application of Hamiltonian and Lagrangian methods 
in fusion plasma physics lies with the Hamiltonian dynamics of magnetic field lines and the 
existence of closed magnetic surfaces [ID] . 

Next, Morrison reviewed the variational structure of the exact Vlasov equation, begin- 
ning with the pioneering work of Low JTJ and culminating with recent works by Cendra 
et al. [T2] and Brizard [13] . and the Hamiltonian and Lagrangian formulations of several 
exact and reduced fluid models (e.g., ideal MHD with flows [H] and reduced MHD [To]). 
Morrison concluded his presentation with a functional-integral derivation of the fluctuation 
spectrum of a Vlasov-Poisson plasma, again highlighting the usefulness of the Hamiltonian 
formalism. 

Allan Kaufman reviewed Boghosian's work JDJ on the covariant Lagrangian formula- 
tions of guiding-center and oscillation-center plasmas. A covariant variational principle is 
formulated for Vlasov-Maxwell plasma dynamics in terms of electromagnetic fields/potentials 
in four-dimensional space-time (with coordinates r M ) and particle orbits in eight-dimensional 
phase space (with coordinates z a = (r^,^)). The advantages of a covariant formulation 
include a frame-independent representation in which polarization and magnetization effects 
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are treated in a unified manner. The covariant Hamiltonian and Lagrangian formulations 
of charged-particle dynamics proceed as follows. First, the covariant Hamiltonian struc- 
ture of charged particle motion in electromagnetic fields is given in terms of the covariant 
Hamiltonian H(z) = p' 1 p ll /2m = —mc 2 /2 and the covariant Poisson bracket 

{zoo, 9 (z)} = d 4^ d 4^- - + - d 4^F, u {r) dg{z) 



dr^ dp^ dr^ c dp^ dp u 

where f(z) and g(z) are arbitrary extended-phase-space functions, and F^ v {r) = d fM A u (r) — 
d u A^(r) denotes the electromagnetic field. We note here that the electromagnetic field 
appears only in the Poisson bracket (|1J), which is derived from the covariant phase-space 
Lagrangian V = (p^ + e A^/c) dr^, and that the covariant Hamiltonian is actually a Lorentz 
invariant |17j . 

The covariant perturbation analysis now proceeds as follows. First, for slowly- varying 
fields i^(ow (with Em) ■ B( ) = 0), Lie-transform methods lead to a covariant guiding- center 
(gc) description (r — > T gc and H — > H gc ). Next, a small-amplitude eikonal wave 

= F {0)flu + e [i{k tl A v -k v A^ exp^e" 1 ^ + c.c. ] + 0(e 2 ) 

is introduced, where e 1 is the eikonal small parameter while the eikonal amplitude A^, 
the eikonal phase 9, and the wave four-vector = e~ 1 d fJ- 8 are slowly-varying functions of 
the space-time guiding-center coordinates, so that the covariant guiding-center phase-space 
Lagrangian becomes T gc — > T gc + eT w . The wave fast time scale is asymptotically elimi- 
nated through a near-identity phase-space transformation from guiding-center coordinates 
to oscillation-center (oc) coordinates, expressed as asymptotic series expansions in powers 
of e and denoted Z(z,e). The phase-space transformation also induces a transformation 
on the Hamiltonian (H gc — > K oc = T~ x H gc = K + e 2 K 2 ), where the covariant quadratic 
ponderomotive Hamiltonian for oscillation-centers is 

K 2 (Z; A,k) = A;(R)K^(Z; k)A u (R). (2) 

Here, the phase-space 4x4 tensor K^ U (Z\ k) appears in the definition of the linear suscep- 
tibility 

X ^(x) = - Jd 8 Z6\R-x) f oc (Z) K^(Z;k), 

where f oc {Z) denotes the oscillation-center Vlasov distribution, which in turn leads to the 
so-called K — \ theorem jTH] 

S XT RIW _ d 2 K 2 
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which relates the oscillation-center ponderomotive Hamiltonian to the linear susceptibility. 
Lastly, Kaufman emphasized the role played by the Noether method in deriving conserva- 
tion laws relating the self-consistent evolution of wave, background plasma, and background 
field [El EH EH- 
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Diego del-Castillo-Negrete presented work on the role of self-consistent chaos in mean- 
field Hamiltonian models of fluids and plasmas. Mean-field Hamiltonian theories are inter- 
mediate between test-particle Hamiltonian theories, in which particles interact with wave 
structures but not with each other, and exact Hamiltonian theories, in which self-consistent 
fields produced by particles are taken into account in wave-particle interactions. Mean-field 
Hamiltonian models can be used to describe the collective dynamics of marginally stable 
fluids (e.g., the dynamics of localized vortices with negative and positive circulation in 
shear flows) and plasmas (e.g., the self-consistent evolution of localized electron holes and 
clumps in phase space |2*T] I22|). During his talk, del-Castillo-Negrete presented the Hamil- 
tonian formulation of the single-wave model (SWM) equations describing the interaction 
of iV particles with phase-space coordinates {(xj, Uj);j = 1,...,N} with a single wave 
(characterized by a time-dependent complex- valued amplitude a). Here, the Hamiltonian 
for the particles is 



N 



k=l 



y - a(t) exp(ix k ) - a*(t) exp(-ix k ) 



(3) 



so that Hamilton's equations dxj/dt = dH/duj and duj/dt = —dH/dxj describe the 
particle dynamics, while the wave dynamics is represented by the equation 

a< ~dl ~ iU a = ~N ^ Ffc exp(_iXfc )' ( 4 ) 

where (a, U, Tf.) are model parameters (e.g., clumps and holes are labeled by > and 
Tfc < 0, respectively). The SWM equations ©-© have their origin in the study of the 
beam-plasma instability [SB]- More recently, the system has been systematically derived 
under more general conditions as a basic model describing the weakly nonlinear dynamics 
of marginally-stable Vlasov-Poisson plasmas and localized vorticity perturbations in shear 
flows [21]. Also, the SWM bears interesting similarities with Hamiltonian models used in 
the study of long-range interacting systems in statistical mechanics [23] ■ 

Numerical simulations based on the single- wave model equations (jHJ), with finite and 
infinite N, show the existence of coherent, rotating dipole states. The coherence of the 
dipole is explained in terms of a parametric resonance betwen the rotation frequency of the 
macroparticles and the oscillation frequency of the self-consistent mean field [22] • The role 
of self-consistent chaos in the formation and destruction of phase-space coherent structures 
was emphasized by del-Castillo-Negrete. Lastly, for some initial conditions, del-Castillo- 
Negrete showed that the mean-field exhibits a self-consistent elliptic-hyperbolic bifurcation 
that leads to the destruction of the dipole and violent mixing of the phase space (see 
Ref. [23] and Section V on Lagrangian Chaos). 



III. LIE-TRANSFORM METHODS 



The Lie-transform methods traditionally used in Hamiltonian perturbation theory in- 
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volve near-identity (or infinitesimal) canonical transformations. Practical applications of 
such perturbative Lie-transform methods normally involve truncation of asymptotic expan- 
sions in powers of a small dimensionless parameter e at finite order and, consequently, 
involve moderate algebraic complexity. Non-perturbative Lie-transform methods, which 
involve finite canonical transformations, are used in Hamiltonian orbit theory (associated 
with accelerator design, for example) and are, therefore, typically associated with tremen- 
dous algebraic complexity (to the point of requiring the use of symbolic manipulation 
algorithms). The present Section focuses on Alex Dragt's presentation on applications of 
non-perturbative Lie-transform methods to the technological challenges associated with the 
construction of powerful particle accelerators |27] , high-resolution electron microscopes 
[2H1, and light optical devices fMW\. 

The Lie-algebraic tools used in accelerator physics and geometric optics include sym- 
plectic matrices, symplectic maps, and Lie transformations. A matrix M is symplectic if it 
satisfies the condition M T • J ■ M = J, where J represents the canonical symplectic matrix; a 
map M. is symplectic if its Jacobian matrix is symplectic; and a (finite) Lie transformation 
is defined as 

(exp:£)$ = 9 + {/, g} + \{f, {/, g} } + (5) 

for any phase-space functions / and g (using Dragt's notation :f:g = {/, g} for Lie opera- 
tors) . 

With these Lie-algebraic tools in hand, we now contemplate their applications. First, 
any symplectic map A4 / can be factorized as 

M f = (exp:f 2 :) (exp:f 3 :) (exp:f 4 :) • ■ ■ , 

where f n is a homogeneous polynomial of degree n in the phase-space coordinates, while its 
inverse map is defined as .Mj 1 = ■ ■ ■ (exp — :f 4 :) (exp — :f 3 :) (exp — :f 2 :). Here, factorization 
is facilitated by the Baker-Campbell-Hausdorff (BCH) theorem 

(exp:f:) (exp:g:) = (exp:h:), (6) 

which states that a new Lie operator :h: can constructed from the Lie operators :f: and :g: 
through the BCH formula 

h = f + g + -{/,<?} + {/, g}} + {g, {g, /}}) + •■■ 

Proof of the BCH theorem (JHJ) relies on the definition Q of the Lie transform and the fact 
that the Poisson bracket { , } satisfies the antisymmetry property and the Jacobi identity. 

Next, the Hamiltonian orbits of charged particles moving in an electromagnetic field 
(or light rays) propagating through an assembly of beam-line elements (or lenses) are 
generated by symplectic maps (i.e., finite canonical transformations). Since the guiding 
action of electromagnetic fields has strong similarities with the guiding action of lenses, the 
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Lie-algebraic tools defined above can thus find applications both in accelerator physics and 
in geometric optics. First, a series of elements (electromagnetic multipoles or lenses) can be 
represented by the map product M. = A4f M. g ■ • -, where each element is represented by 
a symplectic map Aif, with (exp^:) describing paraxial behavior, (exp^:) and (exp-.f^.) 
describing second- and third-order aberrations, etc. Lastly, the product of maps for two 
successive elements can be factorized using the BCH theorem (JUJ). Consequently, a map 
action M. : z m — > z out can be explicitly constructed (by drawing tremendous advantage 
from symbolic manipulation algorithms). 

Applications of Lie-transform methods to the development of electron microscopes 
and linear colliders, for example, can lead to the cancellation of spherical aberration (al- 
ways present in axially symmetric electromagnetic systems) by the introduction of a sex- 
tupole beam-line correcting element. This cancellation can be explicitly calculated by 
Lie-transform methods and a consequent spot-size reduction by at least one order of mag- 
nitude can be observed (e.g., Dragt presented a case where the spot-size was reduced from 
about 3 A to less than 0.1 A). Similar applications in accelerator physics, for example, pro- 
vide accelerator designers with the capability of following charged particles through more 
than one millions of turns around a storage ring. 



Lagrangian (or variational) formulations of fluid and plasma physics provide powerful 
guiding principles for the derivation of exact and reduced dynamical equations and the 
construction of their associated invariants and conservation laws. The present Section 
summarizes the talks presented on applications of Lagrangian methods in fluid and plasma 
physics. 

Before proceeding with our summary, we point out that near-identity transformations 
also exist in fluid dynamics (e.g., ideal magnetohydrodynamics). Here, near-identity trans- 
formations are expressed in terms of the Lie operator T e defined as 



where x( x ; t) is an arbitrary fluid scalar field and the vector field £ represents the ideal 
fluid displacement from a reference (or unperturbed) position x to the exact (or perturbed) 
position x + e£(x, t) of a fluid element. Since near-identity transformations are invertible, 
we find the inverse operator T~ l = exp(— e£- V). For Hamiltonian fluid dynamics (with 
unperturbed fluid velocity u ), the unperturbed evolution operator is do/dt = d/dt + u • V 
while the perturbed evolution operator is d e /dt = d/dt + u e • V, where the perturbed fluid 
velocity can be written as 



IV. LAGRANGIAN METHODS 



exp(e £• V) x(x,t) 



(7) 
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Additional expressions for perturbed fluid quantities, which can also be derived by Lie- 
transform methods, include the perturbed fluid density 



p e (x,t;£) = T7Vo(x,t) detfv^x 



) 



(9) 



the perturbed fluid entropy 

s e (x,t; C) = T7 1 s (x,t), 
and the perturbed (magnetic) vector potential 



(10) 



A e (x,t;£) = T^Aofot) - V (T £ - 1 A (x, t)) -5 e x, 



(11) 



where (5 e x = T e : x — x and we have omitted the gauge term V(T e 1 A • 6 t x) in Eq. (JUJ. 
The expressions i p can also be derived from the Lie derivative of appropriate 

differential forms, e.g., p t d 3 x = exp(— e Cg) [po d 3 x}. Note that the constrained variations of 
the fluid and electromagnetic fields 5ip a = (5p, 5u, 5s, 5 A) used in the standard variational 
formulation of ideal magnetohydrodynamics (HI EH] are defined from ip^ = (p e , u e , s e , A e ), 
as given in Eqs. (|H|)- (fTT|) . as 8i/; a (£) = (dip^/de)\ e=0 . For example, (du e /de) e=0 = (d t + 
u o ■ V)£ — $, • Vuo = 5u yields the well-known expression for the constrained (Eulerian) 
variation of the fluid velocity u |31j . 

Darryl Holm presented a variational principle for Euler-Poincare equations of geophys- 
ical fluid dynamics amenable to asymptotic analysis based on the procedure of Lagrangian 
averaging (221- Holm pointed out that the challenge of geophysical fluid modeling is to 
model large-scale (or slow) dynamics while retaining the effects of small-scale (or fast) 
dynamics, as represented by fluid turbulence. This situation is reminiscent of the single- 
particle case where fast time scales (e.g., orbital time scales associated with gyromotion 
or bounce motion of a charged particle in a strong magnetic field) are asymptotically 
eliminated by Lie-transform methods pQ. In the present case, two averaging procedures 
(Eulerian and Lagrangian) offer opposing characteristics. Here, the Lagrangian average of 
a fluid quantity x(x, t), denoted as x( x , t) and taken at constant Lagrangian coordinate x , 
is defined as x( x , t) = (x(x+ e£, t)), where ( ) denotes an Eulerian average (performed 
at a fixed spatial location x) and x + e£(x, t) represents the exact position of a fluid ele- 
ment whose mean position is x (we note that our notation is slightly different than that of 
Ref. |S2])- Hence, Eulerian averaging (taken at a fixed spatial location) commutes with the 
spatial gradient and partial time derivative and, consequently, the momentum-conservation 
form of hydrodynamical equations is invariant under Eulerian averaging. Lagrangian aver- 
aging (or path- aver aging taken by following a fluid element) commutes with the advective 
time derivative moving with the fluid and, consequently, Lagrangian averaging preserves im- 
portant conservation laws (e.g., Kelvin circulation and potential vorticity). Unfortunately, 
Eulerian averaging does not preserve these conservation laws while Lagrangian averaging 
does not commute with spatial gradients and is (by definition) history dependent. Never- 
theless, a generalized Lagrangian mean (GLM) theory for Euler-Poincare equations can be 
developed within the context of constrained variational principles |H2] - 
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The description of large-scale (slow) dynamics while retaining the effects of the small- 
scale (fast) dynamics is provided by a wave, mean flow interaction (WMFI) model, con- 
structed by asymptotic analysis based on two dimensionless parameters (a, e) associated 
with wave amplitude (a) and eikonal representation (e) of the wave fields. For this pur- 
pose, we replace e£(x, t) in Eqs. ([7J)-(fT0)) with a£(ex, et) exp[ze _1 ^(ex, et)] + c.c. so that 
the perturbed fluid fields p e , u e , and s e are replaced with P( a ,e), u (a,e), and S( a , e )> respec- 
tively. These expressions are then inserted into the variational principle for geophysical 
fluid dynamics and expanded in powers of a and e. The mean-flow Lagrangian is obtained 
at zeroth order in a and e, while at order a 2 , terms representing ray-optics, self-modulation 
(dispersion); WKB stability can be collected at zeroth order, first order, and second order 
in e, respectively. Some approximations of the WMFI variational principle include jS2] : (i) 
the derivation of the mean-flow dynamics by setting a = and varying with respect to 
the mean fields; (ii) variation with respect to the wave fields only yields equations suitable 
for linearized spectral analysis (with unperturbed mean- flow dynamics included); and (iii) 
variation with respect to mean fields and wave fields, after phase averaging is performed, 
yields the WMFI dynamics. 

Lazar Friedland presented his work on an averaged variational principle for multi-phase 
nonlinear waves excited and controlled by synchronized nonlinear mode conversion by slow 
(adiabatic) passage through resonance [331 • The analysis focused on small periodic pertur- 
bations of periodic solutions of the Sine-Gordon (SG) equation jM] u tt (x,t) — u xx (x, t) + 
sinu(x,t) = ef(x,t). By using an eikonal representation for the perturbation, Friedland 
showed that, when the amplitude of the driving perturbation is above a certain threshold, 
one can excite a single phase nonlinear waveform synchronized with the drive. Furthermore, 
as the driving parameters vary in time and/or space, the emerging phase-locked waveform 
slowly evolves near the solution space of the unperturbed problem and, if desired, becomes 
highly nonlinear. The oscillatory solutions of the nonlinear SG equation are the plasma and 
breather oscillations and Friedland showed how each nontrivial nonlinear solution of the SG 
equation created from zero and controlled by a small forcing. Lastly, Friedland presented 
an outline of the theory and numerical simulation results of this excitation process based 
on multi-phase Whitham's averaged variational principle. 

Eliezer Hameiri presented a review of the variational formulation of the Energy- Casimir 
stability analysis of plasma equilibria with general flows j^S]. Beginning with a single- 
fluid model (Hall-MHD), Hameiri identified two stream functions ip an d v 9 associated with 
axisymmetric plasma equilibrium. The first stream function if) is the familiar magnetic 
flux function defined from the magnetic equation B ■ Vif) = while the second stream 
function (p defines flow surfaces through the relation u • Vy? = 0. We note that a general 
plasma equilibrium with flows requires that its entropy s satisfy u ■ Vs = 0, i.e., the 
entropy s((p) is a function of the flow stream function (p. Moreover, Hameiri showed that 
constant-entropy surfaces are also defined from an application of a modified version of 
Ertel's theorem, which leads to the condition B* ■ Vs = 0, where B* = B + a -1 V X u is 
related to the curl of the canonical momentum of the single- fluid (with a = ei/rrii). Hence, 
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in an axisymmetric plasma equilibrium with flows, we may write B = V£ X + / V£ 
and B* = VC X Vy? + I* V£, where I and I* are functions of ip and ip (our notation). 
After identifying all the constants of the motion for a given plasma equilibrium with flows, 
a variational principle is developed jHHl; which yields a description of plasma equilibrium 
with flows in terms of two coupled Grad-Shafranov equations for ip and (p. Hameiri also 
discussed extensions of this variational formulation of equilibria with flows to the multi- 
species magnetofluid case and the classical fluid case. 

Jean-Luc Thiffeault presented work on a unified variational formulation of the nonlin- 
ear stability analysis of static jSZj and stationary |38J plasma equilibria. In particular, 
the methods of Eulerianized Lagrangian Displacements (ELD) [3T| [see Eqs. (J8|)-(|ll|)]. Dy- 
namical Accessibility (DA) j^J, and Energy- Casimir (EC) |H| are shown to be equivalent 
in assessing the stability of magnetohydrodynamic (MHD) equilibria. Thiffeault showed 
how sufficient conditions for the stability of MHD equilibria can be obtained by using the 
Dynamical Accessibility method [4 j. The use of dynamical accessibility ensures that the 
physical perturbations preserve the natural constraints of the system imposed by the Hamil- 
tonian structure. Here, the first and second variations of the fluid fields ( a are expressed 
as S( a = {£?, C a } and 5 2 ( a = \{Q, {(?, ( a } }, where Q = fC a Xa d 3 x is a generat- 
ing functional and { , } is an appropriate Poisson bracket structure. From the ansatz 
£ = P Vxi — X2 Vs + B x V X x 3 , Thiffeault showed that the ELD and DA methods are 
equivalent if the functions xi> X2, and X3 span the space of £ (and vice versa). Lastly, 
Thiffeault commented on the fact that this procedure is a generalization of the work of 
Newcomb (31] and Arnold (2^1 , among others, who make use of Lagrangian displacements. 

V. LAGRANGIAN CHAOS 

Even in the presence of dissipation, a knowledge of Hamiltonian dynamics improves our 
understanding of the role played by dissipative dynamics. John Finn presented numerical 
results of work on self-consistent Lagrangian chaos (i.e., chaotic advection) in locking bi- 
furcations in two-dimensional shear flows. Here, two-dimensional shear flows are modeled 
by the velocity field v = V0 X z and the vorticity scalar field w = ?'Vxv = - V 2 in 
terms of the scalar field <p(x, y, t). Two-dimensional fluid motion x = v is thus described as 
a Hamiltonian system x = dcj)/dy and y = — d<p/dx, with acting as the Hamiltonian. As 
is well known, the fluid motion is integrable if the Hamiltonian is time independent, or can 
be made time independent by a change of variable, e.g. a change of frame; if the Hamilto- 
nian is time dependent, Lagrangian chaos is possible. Lagrangian chaos is defined by the 
property of fluid elements following chaotic trajectories while Eulerian chaos is defined in 
terms of turbulence at a fixed location; as a result of Lagrangian chaos, if a passive scalar 
field is injected into a chaotic flow, the passive scalar field is advected, stretched, and folded 
by the chaotic nature of the flow jlU]- Finn pointed out that Lagrangian chaos does not 
imply Eulerian chaos (although it enhances transport) and that even laminar periodic flows 
can lead to complex Lagrangian trajectories with statistical properties similar to turbulent 
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flows. 



Next, Finn investigated the self-consistent Lagrangian chaos in a fluid system consisting 
of a nonlinear wave due to a shear-flow (Kelvin-Helmholtz) instability and an externally 
imposed sinusoidal perturbation; the fluid was taken to be two-dimensional and incom- 
pressible. Self-consistent Lagrangian chaos was observed in the sense that the vorticity 
was actively advected with the flow jHj. The Navier-Stokes system describing the fluid 
was studied by direct numerical simulation (involving symplectic integration techniques). 
In the region of parameter space studied (moderately large Reynolds number), Lagrangian 
chaos in the presence of smooth flows is relevant, and in this range the bifurcations which 
occur are low dimensional. The main features observed during direct numerical simulations 
are locking and unlocking bifurcations, the latter related to the self-consistent Lagrangian 
chaos j32]- Lastly, Finn discussed the relation with self-consistent Lagrangian chaos in 
shear flows supporting two Kelvin-Helmholtz modes of very different phase velocities [IB] 
(where locking is impossible). 



VI. SUMMARY 



The Mini-Conference was successful in bringing together experts from a variety of fields 
in physics who actively apply Lagrangian and Hamiltonian methods in their research. In his 
talk, Morrison highlighted the surprisingly prevalent lexicon of finite-dimensional Hamilto- 
nian systems (that include such terms as integrability, Poincare sections, KAM invariant 
tori, intrinsic stochasticity, and Chirikov's overlap criterion). Indeed, such concepts are 
now part of our language, independently of whether we use Lagrangian or Hamiltonian 
methods in our research. Kaufman also reminded us of the story of the discovery of the 
third (drift) adiabatic invariant by Northrop. According to Kaufman, Edward Teller called 
in the Sherwood theory group, and posed the problem of how only 2 (known) invariants 
were insufficient to explain the confinement of individual particles in the recently discovered 
van Allen belt. As Kaufman recalled, Ted Northrop on his own deduced the drift invariant, 
and then co-authored a paper with Teller pHj. 

The Mini- Conference also included a poster session during which many presentations 
were made on applications of Hamiltonian and Lagrangian methods in fluid and plasma 
physics. Presenters from France, Sweden, Germany, the United Kingdom, Russia, and 
the United States provided ample evidence of the vitality of this field of research in fluid 
dynamics, laser-plasma interactions, space plasma physics, and theoretical fluid and plasma 
physics. 
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